Critical values for a non-attractive lattice gas model 
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Abstract 

We investigate numerically the critical behaviour of a one- 
dimensional non-attractive lattice gas model that is the continuous- 
time version of the Domany-Kinzel cellular automaton in one of its 
parameter subspaces. The model shows a phase transition in the 
directed percolation universality class of critical behaviour, as ex- 
pected. 

Keywords: critical values, lattice gas, finite-size scaling, Domany- 
Kinzel cellular automaton, directed percolation 



PACS numbers: 02.50.Ga, 05.40.+j, 64.60.Ht 



The Domany-Kinzel (DK) probabilistic cellular automaton (PCA) [|T| is one of the 
most studied PCA in the physics literature, because it is the most general left-right 
symmetric one-dimensional PCA, and has the interesting property of having the mixed 
site-bond directed percolation (DP) process on the square lattice as one of its instances. 
Its defining rules are given in table |l|. The mixed site-bond DP process is given by 
assigning a probability s G [0,1] for a site to be present in the lattice, and a probability 
b G [0, 1] for a bond to exist between any two sites on the lattice. The mixed DP problem 
consists in finding the values of s and b for which an infinite cluster of sites connected by 
the bonds occurs such that one can walk unidirectionally in it, say to the south and to 
the east indefinitely. In the DK PCA this problem is obtained by choosing the transition 
probabilities x — 0, y — sb, and z = sb{2 — b). For s = 1 one obtains the pure bond DP 
problem, whilst for b = 1 one obtains the pure site DP problem. 

In this paper we investigate numerically the critical behaviour of a continuous-time 
one-dimensional non-attractive lattice gas for which some lower bounds on the critical 
points of the PCA version was given recently 0. The model is related with the DK PCA 
in one of its parameter subspaces, and although the model is not on the mixed site-bond 
DP parameter subspace of the DK PCA, it presents DP exponents, as expected on the 
basis of the DP conjecture ||. 

Let rii(t) G {0,1} denote the occupation number of the site £ G A at the integer 
instant t, with AcZa finite lattice with |A| = L sites and periodic boundary conditions 
£ + L = i. The model we are interested in is the continuous-time version of the PCA 
defined by the rules 

n (t + 1) = I (^-i^) + mod 2 witn probability p, 

1 with probability 1 — p. 

The rule table for this PCA is given in table |2|. From tables [l] and |2] we see that our PCA 
is equivalent to the DK PCA with rates x — 0, y — p, and z = 0. We thus see that unless 
we take the unphysical value b = 2 in the site-bond DP subspace of the DK PCA, this 
model does not belong to that subspace. 

Our approach in constructing the continuous-time version for the above PCA is to 
take its non-diagonal transitions, i.e., those transitions for which the final state differs 
from the initial state, and associate with them a stochastic lattice gas with transition 



Table 1: Rule table for the DK PCA. The first line gives the initial neighbourhood, the 
other two lines give the probability at which the state listed at left is reached by the 
central bit. 
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Table 2: Same as table [T] for the PCA defined by Eq. (]T|). 
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rates given by the original PCA rules. This approach has been used before in the PCA 
literature ||, and is equivalent to the so-called 'Hamiltonian' or 'strong anisotropic' limit 
for the transfer matrixes of equilibrium lattice models ||. 

As is well known || , we may write the master equation for interacting lattice g ases as 
a Schrodinger-like equation in Euclidean time, 



dt 



\P(t)) = -H\P(t)), 



(2) 



with \P{t)) the generating vector of the probabilities P(n,t) = (n\P(t)) of observing 
the configuration n = (rii,ri2, ■ ■ ■ , ^l) E {0, 1} A at instant t, and with the infinitesimal 
generator H of the Markov semigroup playing the role of the Hamiltonian. For the non- 
diagonal transitions of the PCA defined by Eq. (|IJ), the operator H can be written as 



H 



(3) 



with the three-body stochastic transition matrix Ht-m+i given by 
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where the three-site basis vectors are ordered as usual, (0,0,0) -< (0,0,1) -< ••• -< 
(1,1,0) -< (1,1,1), and the dots indicate null entries. Proper tensorization of the above 
three-body matrix with unit matrixes in order to obtain the full matrix H is understood. 
Notice that the elements in the columns of i^-i^^+i (and consequently of H) add to zero 
due to the conservation of probabilities, and that its non-diagonal elements are positive, 
since < p < 1. Identifying a particle with the up spin state and a hole with the down 
spin state in the a z basis, the transition matrix H^_i^ + i above is seen to be equivalent 
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to the non-Hermitian quantum spin operator 



He-iAW = - !) I 1 + (! - V)o\ + po^of o£ +1 ] > ( 5 ) 

where a x and a z are the usual Pauli spin-| matrices. The transition matrix (^) resembles 
the analogous matrix for the basic contact process, but with non-standard rates and with 
the elementary process 101 — ^ 111 lacking. This lack is the root of the non-attractiveness 
of the process. (Loosely speaking, attractive interacting particle systems present a ten- 
dency for clustering, as it occurs in ferromagnetic models or in the basic contact process. 
The precise mathematical statement of attractiveness can be found in 0.) 

The lowest gap in the spectrum of H may be used to perform a finite-size scaling 
analysis in the same way as one does in equilibrium problems ||. Around the critical 
point p^p* , the correlation lengths of the infinite system behave like 

^KiiK{p-p*)-»\\K{p- P *)-^, (6) 

where £11 and £j_ are the correlation lengths respectively in the time and space directions, 
v\\ and u± are the corresponding critical exponents, and z = v\\/v± is the dynamical critical 
exponent. For finite systems of size L we expect that 

q^L-^dp-pUL 1 ^), (7) 

where zl and u± } l are the finite versions of z and u±, and is a scaling function with 
$(m ^> 1) ~ On general grounds one expects lim^oo p* L , zl, v± ) l = p*,z,u±. From 
Eqs. (H) and (|7|) we obtain 

^ [iuipD/iwApD] = in [i\\Mpl)/i\Upl)\ = (R) 

ln(L/L') ln(L"/L) Zh " K) 

which through the comparison of three different system sizes V < L < L" furnishes 
simultaneously p* L and zl- Of course, and the gap - Ef ] = of H are related 

We calculated the gaps of H with the power method, which requires only matrix-by- 
vector multiplications that can be carried out efficiently and does not require a diago- 
nalization in the usual, 'QR' sense, a step that may lessen the quality of the data. The 
version of the power method we use takes advantage of the presence of absorbing states, 
and is also suitable for the investigation of time dependent properties of Markov chains || . 

Our results for p* and z are summarized in figure [I]. Curiously enough, despite the 
translational invariance of the lattice gas rules the finite-size estimates behaved better 
for triplets of lengths of the form L',L,L" = 21 - 1,2/, 2/ + 1, I e N. Both sets of 
data behaved irregularly with the system sizes, preventing us from applying the usual 



extrapolation algorithms [10, 11] to them. We are presumably in the presence of strong 
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Figure 1: Finite-size data for the critical point p* and the dynamical critical exponent z 
of the model defined in Eq. (Q). The dashed line in the graph for p* L is the least-squares 
linear fit to the data, whilst the dashed line in the graph for z L represent the best known 
value of zdp- 



finite-size effects and corrections to scaling. Notice that the last three points of the data 
seem to be converging monotonically, but unfortunately we were not able to go beyond 
L = 22 in our diagonalizations. The L = oo values for p* L were obtained through a least- 
squares fit to the curve xl = x^ + aL' 1 , since the data scale well with L~ l (although it 
does so with L~ 2 also, but with a smaller correlation coefficient), whilst for Zl we only 
estimated the mean value of our data. As expected on the basis of the DP conjecture, 
namely, that the phase transition about a single absorbing state in single-component 
systems with a scalar order parameter and in the absence of internal symmetries should 
be in the DP universality class of critical behaviour |§, our lattice gas shows a DP- 
compatible exponent z = 1.58 ± 0.04. The most precise value of z for the DP universality 



class to date is given by z = 1.580 745 ± 0.000 010 ||12|| . The critical point of the model 
is estimated as p* = 0.926 ± 0.004 (LS correlation coefficient 7 = —0.944). This value 
is slightly higher than the critical value p^ K = 0.82 ± 0.01 of the corresponding point 
(x — 0, y — p, z — 0) in the DK PCA []]]. This shift in the critical point for the lattice 
gas version of the PCA was observed before in a study similar to the present one, where 
the properties of the lattice gas on the line x = 0, y = z (corresponding to the pure site 
DP problem) was investigated 0], and is probably a general feature, since asynchronous 
dynamics tend to be more noisy than synchronous dynamics. 

In summary, we conducted numerical diagonalizations of the infinitesimal generator 
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of the continuous-time version of a non-attractive probabilistic cellular automaton (PCA) 
that is an instance of the Domany-Kinzel (DK) PCA. Although our PCA is not in the 
site-bond directed percolation (DP) subspace of the DK PCA, its continuous-time version 
shows DP critical behaviour. Our finite-size data showed an irregular approach to the 
infinite system limit, and this prevented us from obtaining good estimates of the critical 
values. It would be interesting to study this lattice gas by time-dependent Monte Carlo 
methods in order to obtain more accurate critical values. 
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Paulo (FAPESP), Brazil. 
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